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For two types of quintessence models having thawing and tracking properties, there exist analytic 
solutions for the dark energy equation of state w expressed in terms of several free parameters. We 
put observational bounds on the parameters in such scenarios by using the recent data of Supernovae 
type la (SN la), Cosmic Microwave Background (CMB), and Baryon Acoustic Oscillations (BAO). 
The observational constraints are quite different depending on whether or not the recent BAO data 
from BOSS are taken into account. With the BOSS data the upper bounds of today's values of w 
(= wq) in thawing models is very close to —1, whereas without this data the values of wq away from 
— 1 can be still allowed. The tracker equation of state W(o) during the matter era is constrained to 
be W(o) < —0.949 at 95 % confidence level (CL) even without the BOSS data, so that the tracker 
models with w away from —1 are severely disfavored. We also study observational constraints on 
scaling models in which w starts to evolve from in the deep matter era and show that the transition 
to the equation of state close to «j = — 1 needs to occur at an early cosmological epoch. 



I. INTRODUCTION 

Independent observational data such as SN la P, Q , 
CMB [1, g, and BAO [1, i] suggest that about 70% 
of the energy density today consists of dark energy re- 
sponsible for cosmic acceleration. For the constant dark 
energy equation of state w the recent joint data analy- 
sis based on SN la, CMB, BAO, and the Hubble con- 
stant measurement shows that w is constrained to be 
w = -l.OlStJ!;^?^ at 68 %CL 0. If we use the time- 
dependent parametrization w(a) = + 'Wa{l — a), 
where a is the scale factor normalized as a = 1 today, 
the two parameters wq and Wa are constrained to be 
wo = -IM&t^'o^^l and Wa = 0.14i|5 ?° 0. 

One of the simplest candidates of dark energy is the 
cosmological constant characterized by the equation of 
state w = —1, which is consistent with the current ob- 
servational data. However, if the cosmological constant 
originates from the vacuum energy associated with par- 
ticle physics, there is a huge gap between the theoret- 
ical and observed values [8|. Instead, alternative dark 
energy models with dynamically changing w - such as 
quintessence [lol and k-essence [ll| - have been pro- 
posed (see Refs. [IJl for reviews). 

Quintessence is described by a canonical scalar field 
(j) with a potential 1^(0). In the framework of particle 
physics it is generally difficult to accommodate a very 
light scalar field with a mass of the order of the Hub- 
ble parameter k, lO'^^ eV today [H, |ll|. However 
there have been theoretical attempts to construct viable 
quintessence models in particle physics - especially in 
supersymmetric theories [l^. For example, the Pseudo- 
Nambu-Goldstone-Boson (PNGB) [Hi or axions [iTt 
have the potential of the form V{4)) ~ K^[l ± cos((/)//)] 
with suppressed quantum corrections. 



Caldwell and Linder [18| classified quintessence models 
into two classes, depending on the evolution of w. The 
first class corresponds to thawing models, in which the 
field is nearly frozen by a Hubble friction during the early 
cosmological epoch and it starts to evolve once the field 
mass drops below the Hubble rate H. In this case the 
evolution of w is characterized by the growth from — 1. 
The representative potential of this class is the hilltop 
potential such as V{(j)) = -\- cos(0//)]. 

The second class consists of freezing models, in which 
the evolution of the field gradually slows down because 
of the shallowness of the potential at late times. For the 
inverse power-law potential V{4)) = M^'^'P(I)~p {p > 0) 
there is a so-called tracker solution characterized by 
a nearly constant field equation of state w = — 2/(p -I- 2) 
during the matter era [l9[ . In this case the solutions with 
different initial conditions approach a common trajectory 
(tracker) first, which is followed by the decrease of w 
toward —1. 

In addition to tracking freezing models there is another 
sub-class of freezing models associated with scaling solu- 
tions [l^l . In this case the field equation of state scales as 
that of the background fiuid during most of the matter 
era {w ~ Q). The representative potential of this class is 
1/(0) = Vie"-^i"^/^^p' + T/2e"^2"^/^^p', where Afpi is the re- 
duced Planck mass, Ai and A2 are constants with Ai 3> 1 
and A2 < 1 [2l| . In the early matter era the potential is 
approximated as V{(t)) ~ Vie^'^^'^^^^p', which gives rise to 
the scaling solution characterized by it; = with the field 
density parameter fi^ = 3/Af. At late times the domi- 
nance of the potential V2e~'*'^'^/^^p' leads to the rapid de- 
crease of w relative to the tracker case mentioned above. 
Note that the potential V{(j)) = e-^-^^^^p' [(0 - S)" + A] 
proposed in Ref. [l^ also exhibits a similar property 
to that of the double exponential potential (see also 
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Ref. [2 

In this paper we place observational constraints on 
three types of quintessence models: (i) thawing models, 
(ii) tracking freezing models, and (iii) scaling freezing 
models (see Refs. [23| for related works). This analysis 
covers most of quintessence potentials proposed in litera- 
ture. For thawing models we employ the analytic expres- 
sion of w derived in Ref. [25l.[26| under the approximation 
that 1| < 1 (see also Ref. |27l|). The likelihood anal- 
ysis with the SN la and BAO data was carried out in 
Ref. (see also Ref. [1^). We update the analysis by 
using the latest SN la data (Union 2.1 dataset Q) and by 
adding the data of CMB shift parameters measured by 
WMAP7 We also take into account the recent BAO 
data of the BOSS experiment [sij. 

For tracking freezing models one of the present au- 
thors obtained the approximate analytic formula of w 
expressed in terms of two free parameters [13 (see also 
Ref. ^3)- The likelihood analysis based on the SN la 
and BAO data was performed in Ref. [13 . We show that 
adding the CMB and BOSS data further strengthens the 
constraints on the tracker equation of state w 



(0)- 



For 



scaling freezing models it is difficult to derive an analytic 
expression of w, so we resort to numerical simulations to 
find a viable parameter space. 

This paper is organized as follows. In Sec.|ll]we briefly 
review the procedure to derive approximate analytic ex- 
pressions of w in thawing and tracking freezing models. 
The accuracy of those approximations is also discussed by 
solving the equations of motion numerically. For scaling 
freezing models we show that in some cases it is possible 
to fit the evolution of w by using a specific parametriza- 
tion. In Sec. IIIII we first explain the method of our likeli- 
hood analysis based on the SN la, CMB, and BAO data 
and then we proceed to observational constraints on three 
classes of quintessence models. Sec. lIVl is devoted to con- 
clusions. 



equations of motion are 



4> + 3H<j) + v;^ = , 

Pm + SHp,n = , 



(2) 

(3) 
(4) 



where H = d/ais the Hubble parameter, a dot represents 
a derivative with respect to cosmic time t, p^ = (fP' /2 + 
V{(j)), and = dV/d(j). The pressure of the field is given 
by = 4,y2 - F(0). 

In order to discuss the dynamical system we introduce 
the dimensionless variables 12011 



V 



(5) 



Wc define the density parameters of the field and matter, 
as 



Prf 



2 , 2 

X +y , 



(6) 



which satisfy the relation ft^ + 0„i = 1 from Eq. ^ . The 
field equation of state, w =^ P^/ptj,, can be expressed as 



? 9 

X — y 

^2 _|_ y2 



(7) 



The variables x, y, and A obey the following equations of 
motion M, HI 



-lx{l-x'+y') + ^Xy\ 
-yil + x -y )- —Xxy, 



A' = -\/6a-(r - 1)A^ 



(8) 

(9) 
(10) 



where a prime represents a derivative with respect to 
iV = Ina, and 



II. PARAMETRIZATIONS OF QUINTESSENCE 



r = 



vv. 



(11) 



Quintessence [^ is described by a minimally coupled 
scalar field </> with a potential V{(l)). In addition to the 
field (j) we take into account non-relativistic matter with 
an energy density p,„. The action in such a system is 
given by 



S 



Sri 



-9 ^R-\g^''d^4>d.^-V{^) 

(1) 

where g is the determinant of the metric g^jy, R is the 
Ricci scalar, and Sm is the action for non-relativistic mat- 
ter. In the flat Friedmann-Lcmaitre-Robertson- Walker 
background with the scale factor a(i), the dynamical 



Instead of using the variables x and y, it is also convenient 
to derive the dynamical equations of w and $7^. From 
Eqs. ®-([ini) it follows that 



{1-w) 



-3(1 + w) + X^?,{l + w)n^ 



n' 



(12) 
(13) 
(14) 



Depending on the field potential and the initial condi- 
tions, there are several different cases for the evolution 
of w In the following we shall discuss the follow- 

ing three cases: (i) thawing models, (ii) tracking freezing 
models, and (iii) scaling freezing models. 
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Figure 1. The quintessence equation of state w versus a for the 
potential ^ with (a) //Mpi = 0.5, (jn/f = 0.5 {K = 1.9), 
(b) //Mpi = 0.3, = 0.25 {K = 2.9), and (c) f/Mpi = 0.1, 
(l)i/f = 7.6 X 10"'' {K = 8.2). These cases correspond to 
V^g,^{(f)i) < 0, SO that K > 1. The soUd curves show numerical 
solutions, whereas the bald dashed curves describe the results 
derived from the parametrization (|16|l with f^^o = 0.73. 



A. Thawing models 

For thawing models of quintessence the field (j) is nearly 
frozen in the early matter era because of the Hubble fric- 



tion, so that w is close to —1. A typical example for 
thawing models is the potential V{4>) = A^[l + cos((/)//)]. 
One can regard w = — 1 as the fixed point of Eq. (|12p . 
For X such a point is not stable, and hence w starts 
to grow at the late cosmological epoch. 

If we assume that A is nearly constant, one can express 
w in terms of fi^ by using Eos, and under the 
approximation |1 + lyj <C 1 [23|. This neglects the ef- 
fect of the field mass squared V^^, but it is possible to 
derive a more elaborate form of w with the mass term 
taken into account [2^ [2^. In doing so, the potential is 
expanded around the initial field value </),; up to second 
order, i.e., Viq^) = ELo (</> " Intro- 

ducing a rescaled field u = (cj) — (j>i)o?'/'^ and using the 
approximation ~ —p^ ~ —V{4)i), Eq. ([3]) yields 



4 Ml, 



^^'''V^{cj)^) . (15) 



Provided that |u; -f 1| <C 1, the evolution of the scale 
factor can be approximated by that of the ACDM model, 
i.e., a{t) = [(1 - sm}a^^{t/tA), where U = 

2M.p\/ ^J'SV{(|>i) and VI^q is the today's density parameter 
of quintessence. On using this solution, Eq. ((T5|) is inte- 
grated to give the field in terms of the function of t and 
hence the field equation of state ~ —1 -I- 4)^ /V{(j)i) 
is known. This process leads to the following analytic 
expression of w [2^ [2^ : 



u;(a) =-l + (l + u;o)a3(^-i) 



[K - F[a))[F{a) + 1)^ + (A' + F{a)){F{a) - 1) 



K 



(K - + 1)^ + (a + n-^' 



"4)0 



00 



(16) 



where wq is the value of w today, and 




3F(0,) 



1 + (fi-^i - l)a- 



(17) 
(18) 



The solution ((T5)) is valid for > 0. The equation of 
state (jl6p is expressed in terms of the three parameters 
Wo, fi^o, and K. 

As a concrete example, let us consider the hilltop po- 
tential 



y(</)) = A* [1 -f- cos(0//)] 



(19) 



In this case the parameter (jl7p is given by 

-I 1/2 



1 



Mpi 
/ 



COs((/)i//) 



l + cos(0,//) 



(20) 



If < (j)i/f < tt/2 and 7r/2 < (t>i/f < n, one has K > 1 
and AT < 1 respectively. In the former case the poten- 
tial is approximately given by F(0) ~ 2A''[1 — (fP /{Af'^)] 
around (/) = 0, whereas in the latter case it is approxi- 
mated as V{(t>) sa A4((/) - tt/) V(2/^) around (/> = tt/. 

In Fig. [1] we plot the numerical evolution of w ver- 
sus a for AT > 1 with several different values of / and 

. The bald dashed curves correspond to the results de- 
rived by the analytic expression ([T5|) . which show good 
agreement with the numerically integrated solutions for 
Wq < —0.3. For AT larger than 10, the initial displace- 
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Figure 2. The quintessence equation of state w versus a 
corresponding to the potential (|19p with f /Mp\ = 2. Each 
solid curve represents the numerical evolution of w for (a) 
(t>,/f = 1.8 {K = 0.95), (b) 0,// = 2.1 (A- = 0.81), (c) 
ct>i/f = 2.3 (A' = 0.58), and (d) 0,// = 2.37 (AT = 0.40). In 
these cases V^^^^cjii) > 0, so that K < 1. The dashed curves 
are derived from the parametrization p6|) with Jl^o = 0.73. 



ment of the field is required to be close to to avoid its 
rapid roll down. In such cases the field mass is largely 
negative, which leads to the tachyonic instability of field 
perturbations. If the field reaches the potential minimum 
by today and it starts to oscillate, numerical simulations 
show that Eq. (|16p is no longer reliable. We set the prior 
if < 10 in the likelihood analysis of Sec. IIII Al 

Figure illustrates the variation of w for K < 1 and 
//Mpi = 2 with several different values of (pi. While the 
analytic estimation ([T6)) is trustable for K close to 1, it 
starts to lose accuracy for smaller K. This reflects the 
fact that the field is initially located away from the po- 
tential maximum. For smaller K the field mass becomes 
larger, so that the Taylor expansion around cj) = (pi tends 
to be more inaccurate because of the rapid variation of 
the field. Figure [5] shows that the analytic is reliable 
for 0.5 < if < 1 and wq < -0.8. 



B. Tracking freezing models 

In Eq. (|12p there is another fixed point given by 
3(1 + w) 



A2 



(21) 



along which w is constant. This corresponds to the 
tracker that attracts the solutions with different initial 
conditions to a common trajectory. A typical example 
having a tracker solution is the inverse power-law poten- 
tial V{(f>) = M^+P(j)-P {p > OV The condition under 
which the tracking occurs is jiolpl 



r > 1, 



(22) 



In this case the variable A approaches 0. 

From dH]) it follows that n'^/n^ = -2\' /\. Using this 
relation with Eqs. (fT3|) and (|14|) in the regime il^ ^ 1, 
the field equation of state along the tracker is given by 



W = ^(o) 



2(r-i) 
2r- 1 



(23) 



For the potential V^cp) = M'^+I^-p one has F = 1 + 1/p 
and hence W(o) = -2/(j) + 2) 

The result ([^5]) was derived by neglecting the contri- 
bution of fi^, but its effect can be accommodated by 
dealing with $7^ as a perturbation to the 0-th order so- 
lution ([23|) [32| . For the models in which F is constant, 
the first-order perturbation around Wj-g) obeys 



2d\ 
I -r^OW 
da'' 

9 



■•^^)n4a)=0, (24) 
where ^ltj,{a) is approximated by the 0-th order solution 



-3w(, 



(25) 



Substituting Eq. ([25]) into Eq. ((24| , we obtain the follow- 
ing solution to w{a) (valid for all orders in il^) [3^ : 



^ A similar condition for the k-essence Lagrangian P{(j>, X) 
V{iP)W{X) with X = 02/2 is given by T > 3/2 H. 
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= -(0) + E ._.Ln..,.„.^ ! o...i° ... ) (26) 



1 - (n + l)w(o) + 2n(7i + l)^^^^ V 1 " ^^(a) . 
(l-^(o)Ko) ^ , , , (l~"'(o)Ko)(8^>(o)-l) 2 



2(1 - wl^^)wf^^i'iw^o) - l)(18«;(o) + 1) 



n^iay + ---, (27) 
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Figure 3. The quintessence equation of state w versus the 
scale factor a for the inverse power-law potential V{(fi) — 
M^ff)-^. The solid curve shows the numerical evolution for 
the tracker solution. The other three curves correspond to the 
analytic solution (|27|l with 1-st, 2-nd, and 3-rd order terms 
taken into account. 

where in the second and third lines we car- 
ried out the expansion {V,tf,{a)/{1 — il^{a)))^ ~ 
f^0(a)"(Em=o^<^(a)'")"- The equation of state 1^ is 
expressed in terms the two parameters wjo) and il^o- 

In Fig.|3]we show the 1-st, 2-nd, and 3-rd order tracker 
solutions of Eq. (gT]) for the potential ^(0) = M^cj)-^. 
The approximated formula ^7} tends to approach the 
full numerical solution by adding higher-order terms of 

C. Scaling freezing models 

The scaling solution [2d| is a special case of a tracker 
along which fl^ in Eq. (PT|) is constant with < fJ^ < 1. 
From Eq. it then follows that w = and hence 

^4, = 3/A^ during the matter era. Since A is constant, 
one has F = 1 from Eq. ([Ti)) . This corresponds to the 
exponential potential V{(j)) = Vbe~'^'^/*^p' , where Vq is 



a constan10. In this case, however, the field equation of 
state is the same as that of the background fluid (w = 0), 
so that the system does not enter the phase of cosmic 
acceleration. 

This problem can be circumvented for the following 
model [21| 

T/(0) = y^e-^i0A-fpi + y^e-^^"^/*^"' , (28) 

where and Vi {i = 1,2) arc constants. If Ai ^ 1 
and A2 < 1 then the solution first enters the scaling 
regime characterized by £7^ = 3(1 + !(;,„)/ Af, where w„i is 
the equation of state of the background fluid. From the 
bound coming from big bang nucleosynthesis [13 there 
is a constraint < 0.045 (95% CL) during the radia- 
tion era (wm = 1/3), which translates into the condition 
Ai > 9.4. The scaling matter era {il^ = 3/Af, w = 0) is 
followed by the dark energy dominated epoch driven by 
the presence of the potential ^26"^='*/^^?'. If < 3, the 
solution approaches another attractor characterized by 
flfP = 1 and w = — l-|-A|/3 [l^]. The cosmic acceleration 
occurs for Aj < 2. 

The onset of the transition from the scaling matter 
era to the dark energy dominated epoch depends on the 
parameters Ai, A2, and V2/V1. Numerically we find that 
the transition redshift is not very sensitive to the choice 
of V2/V1, so we study the case V2 = Vi. In Fig. |4] we 
plot the numerical evolution of w for A2 = with three 
different values of Ai . For larger Ai the transition to the 
dark energy dominated epoch occurs earlier. 

It is possible to accommodate the above variation of 
w analytically by using the paramctrization proposed by 
hinder and Huterer [s^l (see also Refs. (39|): 

where Wp and Wf are asymptotic values of w in the 
past and future respectively, at is the scale factor at the 
transition, and t describes the width of the transition. 



^ For the k-essence Lagrangian P{4>, X) (where X = 
—g^^'^ dfi4)dv(j>/2) the condition for the existence of scaling 
solutions restricts the Lagrangian in the form P = X g{Y), 
where g is an arbitrary function in terms of Y = Xe^'^^^^'pi [36l . 
Quintessence with the exponential potential corresponds to the 
choice g{Y) = 1 - Vo/Y. 
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Figure 4. The quintessence equation of state w versus a for 
the potential ^ with (a) Ai = 10, A2 = 0, (b) Ai = 15, 
A2 = 0, and (c) Ai = 30, A2 = 0. The solid curves show the 
numerical solutions, whereas the dashed curves represent the 
results derived from the parametrization (|29|) with Wp = and 
Wf — —1. Each dashed curve corresponds to (a) at — 0.23, 
r = 0.33, (b) at = 0.17, r = 0.33, and (c) at = 0.11, r = 0.32. 



The scaling solution during the matter era corresponds 
to Wp — 0. For A2 = one has Wf = —1, so that 
the parametrization ([^ reduces to w{a) = — 1 + [1 + 
(a/at)^/"^]"^. Figure |4] shows that the parametrization 
(|29p can fit the numerical evolution of w very well for 
appropriate choices of at and r. For A2 = the transi- 
tion width is around r w 0.33, while at depends on the 
values of Ai. In this case one can carry out the likelihood 
analysis by fixing r = 0.33 and find the constraints on at 
(see Ref. [40| for a related work). 

If A2 7^ 0, then the field equation of state finally ap- 
proaches the value ?« = — l-|-A|/3. Numerically we find 
that w tends to have a minimum for larger A2 before the 
solutions reach the attractor. If Ai ~ 20, for example, 
the minimum appears for A2 ^ 0.3 (sec Fig. [5]). In or- 
der to place observational bounds on A2 in such cases, 
we need to resort to numerical simulations without using 
the parametrization 



III. OBSERVATIONAL CONSTRAINTS 

In this section we place observational constraints on 
the three types of quintessence models separately. We 
use the recent SN la data called Union 2.1 dataset 0, 
shift parameters provided by WMAP7 and the BAO 
distance measured by SDSS7 and by BOSS [U. 

In SN la observations the luminosity distance d^lz) = 
(1 -I- z) H~^{z)dz is measured by the difference (dis- 



Figure 5. The quintessence equation of state w versus a for 
the potential with Ai = 20 and A2 = 0.2, 0.5, 1, y/2. 



tance modulus) of the apparent magnitude m(z) and the 
absolute magnitude M, as 



yu(2;) = ■m{z) - M = 5 log^o [dL{z)/lQ pc] 



(30) 



For the observed distance modulus /iobs(^i) with the er- 
rors cr^,i, the chi square of the SN la measurement is 



given by 



2 _ Mobs(^i) ~ Mth(^i) 

AsNIa — 2-^ 



(31) 



^,2 



where ^th{zi) is the theoretical value of ^{zi) known for 
a given dark energy model. 

The position of the CMB acoustic peaks is related with 
the following two parameters 41 1 



7^ = \ a 



>(0) 



H{z)/Ho 



la=^A^- (32) 



Here is the redshift at the decoupling epoch, d\i^ (z*) 



TZ/[HQ\j rim^] is the comoving angular diameter distance 
to the last scattering surface, and rs(zH.) is the sound 
horizon defined by 



rs{z*) 



dz 



)(0) 



(33) 

where ^^^'^ and ^^'^ are the density parameters of 
baryons and photons, respectively. For the redshift z^ 
we use the fitting formula of Hu and Sugiyama [l^l . The 
chi square associated with the WMAP7 measurement is 



XcMB 



'-CMB^^CMB^CMB : 



(34) 



7 



where X^MB = (^a- 302. 09, 71-1.725,2* -1091. 3), and 
the inverse covariance matrix is given by Q 



2.305 29.698 -1.333 
29.698 6825.27 -113.18 
-1.333 -113.18 3.414 



(35) 



In BAO observations the ratio rBAo(^) = 
fsizd)/ Dv{z) is measured, where rs{zd) is the sound 
horizon at which the baryons are released from the Comp- 
ton drag of photons and Dv{z) is the effective BAO 
distance defined by = [{]^ H-^{z)dzf z / H{z)Y/^ 

0. For the redshift Zd we use the fitting formula of 
Eisenstein and Hu [i^. The chi square of the SDSS7 
measurement is given by 



XbAO,SDSS7 ~ ^BAo'-^BAO-^BAO : 



(36) 



where X^AO ^ (rBAo(0.2)-0.1905, tp 
and the inverse covariance matrix is V 



)(0.35)-0.1097), 



"-^BAO 



30124 -17227 
-17227 86977 



(37) 



We also use the BAO data from the WiggleZ and OdFGS 
surveys, for which Awiggioz(2 = 0.6) = 0.452 ± 0.018 
m, and Amfgs{z = 0.106) = 0.526 ± 0.028 HI, where 
A{z) is defined as Athiz) — Dv{z)\JVljnoHfJ z. Finally, 
we use the latest and most precise data point from the 
BOSS experiment [3l|, for which 

1 Abao(^ = 0.57) 13.67 ± 0.22 . (38) 

The total chi-square from the three datasets is 

= xIn la + XcMB + xIao ' (39) 
where the best-fit corresponds to the lowest value of . 



A. Thawing models 

Let us study observational constraints on thawing 
models given by the equation of state ([T6|) . We first fix 
the value of K and vary the two parameters and fi^o- 
Although quintessence corresponds to the case wq > —1, 
we also allow the regime < — 1 in the likelihood anal- 
ysis. In fact it was shown in Ref. [i^l that the equation 
of state of a phantom scalar field can be accommodated 
by the analytic formula ([T51) . The likelihood results are 
quite different depending on whether the recent BAO 
data from BOSS |3l| are included or not, so we present 
two constraints with/without the BOSS data. 

In Fig. [6] we plot observational bounds on wq and fi^o 
for is: = 1.01. In the absence of the BOSS BAO data the 
dark energy equation of state today is constrained to be 
-1.219 < Wo < -0.930 (95 %CL). Dutta and Scherrer 
(25| showed that even the value Wq = —0.7 is allowed 
from the SN la data alone. Adding the WMAP7 and 
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Figure 6. Icr (red) and 2a (green) observational regions in 
the (wo, f2,^o) plane for K — 1.01 constrained by Union 2.1 
01, WMAP7 [ll, SDSS7 01, and the recent BOSS data [3l[. 
The dot-dashed and dotted curves represent la and 2a con- 
tours res pec tively found by removing the BOSS data from the 
analysis [311 ]. 
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Figure 7. The same as Fig. [6l but for K = 4. 



SDSS7 data gives rise to much tighter bounds on wq- 
This is the case for 17^o as well. We obtained the bound 
0.7078 < n^o < 0.740 (95 %CL), whereas the values 
0.68 < n^o < 0.8 are allowed in Ref. [H. 

We also show the likelihood results with the latest 
BOSS data as red {la) and green {2a) regions in Fig. |6l 
We find that Wq and Sl^o are constrained to be —1.242 < 
wq < -0.995 and 0.705 < fl^o < 0.734, respectively, 
at 95 %CL. Rather surprisingly, the allowed parameter 
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-2 -1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6 



Figure 8. la (red) and 2a (green) observational contours in 
the {wo,Q,,j,o) plane marginalized over fi^o- We set the prior 
0.1 < if < 10. The dot-dashed and dotted curves correspond 
to the 1(7 and 2a constraints without the BOSS data, respec- 
tively. 

space in the regime wq > — 1 is very narrow. In particu- 
lar the ACDM model, which corresponds to = — 1, is 
outside the la observational contour. This comes from 
the fact that the BOSS data ([55)) does not allow a large 
parameter space with w > — 1 at the redshift z = 0.57. 

Figure [7] shows observational bounds on wq and fl^o for 
K = i. Without the BOSS data we obtain the bounds 
-1.408 < Wo < -0.855 and 0.711 < f^^o < 0.740 (95 
%CL), whereas with the BOSS data these parameters 
are constrained to be -1.460 < wq < -0.982 and 0.709 < 
fl^o < 0.734 (95 %CL) respectively. In the former case 
the upper bound of ifo gets larger than that for K = 1.01. 
This reflects the fact that w changes more rapidly at 
late times for larger K. In the presence of the BOSS 
data, however, the allowed parameter space in the regime 
u?o > — 1 is still very tiny. 

We also vary the three parameters wq, ^4,0, and K 
in the likelihood analysis with the prior 0.1 < if < 10 
updating the analysis made in Ref. [2^. With the BOSS 
data taken into account, the best-fit model parameters 
are found to be wq = —1.102, O^o = 0.719, and K = 0.1 
with Xmin = 568.59. Under the prior wq > —1 the best-fit 
parameters are wq = —0.992, fi^o = 0.727, and K = 9.95 
with Y^- = 574.86. In the latter case w varies only 
slightly from —1 around today, whose Xmin larger than 
that for the best-fit ACDM model (xacdm = 573.89 with 
n^o = 0.730lt°;™?^ at 95% C.L.). 

Although the case wq < —1 is plagued by a ghost 
problem, the dynamics of w given phenomenologically 
by Eq. (fTB|) is able to fit the data quite well for wq « 
— 1.2. In fact, with two parameters more than the 
ACDM, according to the Akaike information criterion 



(where we should add twice the number of free parame- 
ters to the original x^) 113, the best-fit corresponds to 
= 568.59 + 2 x 3 = 574.59, whereas in the ACDM 
Xacdm = 573.89 + 2 = 575.89. Therefore the model with 
Wq < —1, even with three parameters, can compete with 
the ACDM. 

In Fig. [8] wc show observational constraints in the 
{wo,K) plane marginalized over fi^o- In the regime 
0.1 < if < 1 the constraints on wq are practically in- 
dependent of K, i.e., -1.212 < Wo < -1.003 (95 %CL) 
with the BOSS data, which should be compared with the 
previous result of Ref. [l^] (slightly updating the result 
in Ref. [H), -1.14 < wq < -0.92 for K <2. In the pres- 
ence of the BOSS data the allowed region shifts toward 
w less than —1, as we see in Fig. [S] 

For K > 1 the lower bound on w gets smaller with in- 
creasing K, whereas the upper bound on w is practically 
unchanged, li K = 9.95, for example, wq is constrained 
to be -2.059 < Wo < -1.014 (95 %CL). For K larger 
than the order of 1 the field equation of state can rapidly 
increase in low redshifts, but such rapid growth of w is 
strongly disfavored from the BOSS data. If we do not 
take into account the BOSS data, the growth of w away 
from Wq = —1 can be still allowed. For K > 10 the an- 
alytic expression is not completely reliable because 
of the rapid rolling down of the field along the potential 
with carefully chosen initial conditions. 



B. Tracking freezing models 

Let us proceed to observational constraints on tracker 
solutions whose equation of state is given by Eq. (|27p . 
Although W(o) is theoretically larger than —1 for 
quintessence, we do not put prior W(o) > — 1 in the ac- 
tual likelihood analysis. In Fig. [HI we show la and 2a 
observational contours in the (i(;(o)j ^</>o) plane. 

Without the BAO data, the tracker equation of state 
is constrained to be -1.188 < W(^o) < -0.949 (95 % 
CL). Meanwhile the analysis of Ref. [1^ based on the 
Constitution SN la and the SDSS BAO data gives the 
bound -1.19 < W(o) < -0.90 (95 % CL). Hence the 
upper bound of W(o) becomes tighter by including the 
WMAP7 data of CMB shift parameters. For the po- 
tential V{(j)) = M^+P(l)-P {p > 0) the 2a constraint 
W(^Qj < —0.949 translates into p < 0.107. 

If we take into account the BOSS data in the analy- 
sis, the best-fit model parameters are found to be W(^q) = 
-1.097, and O^o = 0.717 with Xmin = 568.39. In this 
case the Akaike criterion gives = 568.39 + 2 x 2 = 
572.39, which is smaller than the ACDM value Xacdm ~ 
575.89 with the difference more than 2. The 2a observa- 
tional bounds arc found to be —1.211 < i(;(o) < —0.998 
and 0.701 < f^^o < 0.733 (95 % CL). The upper bound 
of W(o) is very close to —1, which shows that the tracking 
quintessence away from —1 is strongly disfavored from 
the data. 

If we put the prior W(o) > —1 in the analysis with 
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Figure 9. Icr (red) and 2a (green) observational contours on 
tracking freezing models in the (wojfl^o) plane. The dot- 
dashed and dotted curves correspond to the la and 2a con- 
straints without the BOSS data, respectively. 

the BOSS data, we find that the best-fit parameters are 
Wf^o) = -0.998 and fl^o = 0.730 with x'Lin = 574.09. We 
also obtain the upper bound W(o) < —0.964 (95 % CL). 
According to the Akaike criterion the best-fit corresponds 
to Xmin = 578.09 and hence the tracker solution with 
w;(o) > — 1 is not favored over the ACDM. 

C. Scaling freezing models 

For the scaling models characterized by the potential 
([28ll with A2 = we already showed in Sec. Ill CI that 
the evolution of w can be well approximated by Eq. (|29p . 
This parametrization admits an exact solution of H/Hq 
in terms of a and the other free four parameters, which 
is numerically convenient for analyzing the data [38| . In 
Fig. [TU] we plot observational bounds in the (at,r2mo) 
plane derived by using the parametrization (j29p with 
Wp = 0,Wf = —1, and r = 0.33. If the BOSS data are 
taken into account in the analysis, the transition redshift 
is constrained to be at < 0.23 (95%CL). The minimum 
of is found to be at = with ilmo = 0.27, i.e., the 
ACDM limit. The case (a) shown in Fig. 2] corresponds 
to the marginal one in which the model is within the 2a 
observational contour. This means that w needs to ap- 
proach — 1 in an early cosmological epoch {w < —0.8 for 
the redshift z < 2). Even without the BOSS data the up- 
per bound on at is practically unchanged, which refiects 
the fact that w is close to —1 at low redshifts. 

For A2 = we also carry out the likelihood analysis 
without resorting to the parametrization (j29|) . Numeri- 
cally we solve the dynamical equations ((5|)- ([TU1) by tun- 
ing initial conditions to find the evolution which gives 
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at 

Figure 10. la (red) and 2a (green) observational contours on 
at and SlmO for the parametrization (|29p of the scahng sohition 
with Wp — 0, Wf = —1, and r = 0.33. The dot-dashed and 
dotted curves correspond to the la and 2a constraints without 
the BOSS data, respectively. 




0.255 0.26 0.265 0.27 0.275 
^mO 

Figure 11. la (red) and 2a (green) observational contours 
on Ai and QmO for the potential (|28|) with A2 = 0. The 
dot-dashed and dotted curves correspond to the la and 2a 
constraints without the BOSS data, respectively. 



the desired values of Qmo Si,nd flro for the potential 
V{(j)) = T4e-^i'^/*^p' + V2. We put the prior Ai > 9.4 
coming from the constraint of big bang nucleosynthesis. 
As long as T^e"^!*^/*^?' ^ V2 the solutions approach the 
scaling fixed point x = y = •\/6/(2Ai) during the matter 
era, so the initial conditions of x and y are irrelevant to 
the likelihood analysis. In Fig. [TT]we show observational 
bounds on the parameters Ai and flmo- The parameter 
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Figure 12. Icr (red) and 2cr (green) observational contours 
on Ai and A2 for the potential (|28p with QmO = 0.269. The 
dot-dashed and dotted curves correspond to the la and 2a 
constraints without the BOSS data, respectively. 

Ai is constrained to be Ai > 13 (95 % CL), which is con- 
sistent with the results presented in Fig. (TU] For larger 
Ai the transition from w = to w = —1 occurs earlier, 
so that the models are favored from the data. 

For non-zero values of A2 we do not have an analytic 
expression of w, so we solve the background equations 
of motion numerically with the priors Ai > 9.4 and 
A2 > 10^**. We set the latter prior because the A2 = 
case was already discussed above. Varying the three pa- 
rameters Ai, A2, and f2,„0j we find that the best-fit model 
parameters arc Ai = 54.94, A2 = 10^**, and il„io = 0.269 
with Xmin — 574.18. Then the Akaike criterion gives 
= 574.18 2 X 3 = 580.18, which is larger than the 
ACDM value Xacdm = 575.89 with the difference more 
than 4. We also obtain the following bounds 

0.262 < Q,no < 0.276 (68% CL), (40) 

0.256 < n,no < 0.279 (95% CL), (41) 

Ai > 16.3 (68% CL), (42) 

Ai > 11.7 (95% CL), (43) 

A2 < 0.361 (68% CL), (44) 

A2 < 0.539 (95% CL). (45) 

Figure [T^] shows the observational constraints in the 
(Ai, A2) plane for ^Imo = 0.269. If the three parameters 
Ai, A2, r^mo are varied in the likelihood analysis, it is dif- 
ficult to marginalize over flmo in the range < f2mO < 1 
because the solutions are prone to numerical instabilities 
around the tail regions of flmo- Hence we use the fixed 
density parameter flmo = 0.269, which corresponds to 
the best-fit value when the three parameters are varied. 

If A2 > 0.5, the model is excluded at 95 % CL. This is 
associated with the fact that w possesses a minimum for 



larger A2 (see Fig. Therefore, in the context of this 
model, the expansion of the Universe accelerates forever. 
If Ai = 20, for example, the parameter A2 is constrained 
to be A2 < 0.3 (95 %CL). In fact the minimum of w 
appears for A2 > 0.3, which leads to the deviation of w 
from —1 today. In summary, the model is within the 
2(7 observational contour provided that A2 ^ 0.1 and 
Ai > 1. 



IV. CONCLUSIONS 

In this paper we placed observational bounds on three 
types of quintessence models: (i) thawing, (ii) tracker, 
and (iii) scaling models. We used the recent data of 
SN la from Union 2.1, the CMB shift parameters from 
WMAP7, and BAG from SDSS7 and BOSS, by which 
the background cosmic expansion history is constrained 
from the distance measurements. 

In thawing models where the field starts to evolve at 
the late cosmological epoch, the dark energy equation of 
state can be expressed as Eq. p6)) with the three param- 
eters Wo, ^rpo- a-nd K. The parameter K is related to 
the field mass at the initial stage. We put observational 
bounds in the {wo,^rf>o) plane for several different val- 
ues of K (see Figs. [6] and[7l). If we take into account 
the BOSS data ([55)1 in the likelihood analysis, the upper 
bounds on wq are very close to —1 independent of the 
values of K ranging in the region 0.1 < K < 10 (under 
which the analytic formula p6)l is trustable). Without 
the BOSS data the deviation of w away from —1 today 
can be still allowed, as seen in Fig. [51 

In tracking freezing models where w is nearly con- 
stant {w w u'(o)) during the matter era, there is the 
analytic formula (j27p derived by considering a homo- 
geneous perturbation around the tracker. Without the 
BOSS data the tracker equation of state is constrained 
to be -1.188 < W(o) < -0.949 (95 %CL), whose up- 
per bound is tighter than wm\ < —0.90 derived by us- 
ing SN la and SDSS7 data Inclusion of the BOSS 
data gives the upper bounds of u'(o) close to —1. We find 
that without the prior on W(o) the constraint is —1.211 < 
W(o) < —0.998 (95 % CL) and with the quintessence prior 
Wfo) > —1 the upper bound is W(o) < —0.964 (95 % CL). 

For the potential (^5)) with Ai » 1 and A2 <C 1, w 
is close to during the deep matter era because of the 
dominance of the steep potential Vie~''*i"^/*^p' . The field 
equation of state starts to decrease after the potential 
Y^^-X2'P/M^i dominates over Vie'^^'^/'^^^K For larger Ai 
the exit from the scaling regime {w = 0) occurs earlier. 
When A2 = we found that the evolution of w can be 
approximated by the parametrization (|29|) with Wp = 0, 
Wf = —1, and T ~ 0.33, where the transition scale factor 
at depends on Ai . Using this parametrization we derived 
the bound at < 0.23 (95 % CL), which translates into the 
constraint Ai > 13. This is consistent with the bound 
shown in Fig. [TO] derived by solving the background equa- 
tions of motion numerically. 
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For the potential ([28]) with non-zero values of A2, the 
field equation of state tends to have a minimum for larger 
A2. As we see in Fig. [121 the models with A2 > 0.5 is ex- 
cluded at 95 % CL. The parameters A2 < 0.1 and Ai > 17 
are allowed from the data, which shows that the early 
transition from the scaling regime to the regime close to 
u" = — 1 is favored. For the scaling models the observa- 
tional constraints on Ai and A2 are not very sensitive to 
inclusion of the BOSS data, because w evolves toward 
— 1 at the late cosmological epoch. 

In current observations there is no strong evidence that 
quintessence is favored over the ACDM model from sta- 
tistical point of view. If we extend the analysis in the 
regime w < —1, we showed that some of the models 
studied in this paper can compete with the ACDM model 
according to the Akaike information criterion. We note 
that in many modified gravity models such as f{R) grav- 
ity , (extended) Galileons H^j , and Lorentz violating 
theories [1^ it is possible to realize w < —1 without 



having ghosts and instabilities. It remains to see how 
much extent future high-precision observations can con- 
strain quintessence and modified gravity models from the 
background expansion history as well as from the cosmic 
growth of density perturbations. 
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